Evidence of Josephson-coupled superconducting regions at the interfaces of Highly 

Oriented Pyrolytic Graphite 
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Transport properties of a few hundreds of nanometers thick (in the graphene plane direction) 
lamellae of highly oriented pyrolytic graphite (HOPG) have been investigated. Current- Voltage 
characteristics as well as the temperature dependence of the voltage at different fixed input currents 
provide evidence for Josephson-coupled superconducting regions embedded in the internal two- 
dimensional interfaces, reaching zero resistance at low enough temperatures. The overall behavior 
indicates the existence of superconducting regions with critical temperatures above 100 K at the 
internal interfaces of oriented pyrolytic graphite. 

PACS numbers: 74.10.+v,74.45.+c,74.78.Na 



Superconductivity in doped graphite goes back to 1965 
when it was first observed in the potassium intercalated 
graphite C 8 K^. Since then a considerable amount of 
studies reported this phenomenon, reaching critical tem- 
peratures T c ~ 10 K in intercalated graphite @, Q and 
above 30 K - though not percolative - in some HOPG 
samples [1] as well as in doped graphite d-d . Theoreti- 
cal works that deal with superconductivity in graphite as 
well as in graphene have been published in recent years. 
For example, p-type superconductivity has been pre- 
dicted to occur in inhomogeneous regions of the graphite 
structure @ or d— wave high-T c superconductivity fTol| 
based also on resonance valence bonds [111], or at the 
graphite surface region due to a topologically protected 
flat band[l2j. Following a BCS approach in two dimen- 
sions (with anisotropy) critical temperatures T c ~ 60 K 
have been estimated if the density of conduction electrons 
per graphene plane increases to n ~ 10 cm -2 , a den- 
sity that mig ht be induced by defects and/or hydrogen 
ad-atoms[13j or by Li deposition 14]. Further predictions 
for superconductivity in graphene support the premise 
that n > 10 13 cm~ 2 in order to reach T c > 1 KflE l [l6 | . 

The possibility of high-temperature superconductiv- 
ity at surfaces and interfaces has attracted the atten- 
tion of the low-temperature physics community also 
since the 60's[17|. Recently, superconductivity has been 
found at the interfaces between oxide insulators [l8[ as 
well as between metallic and insulating copper oxides 
with T c > 50 K[l!|. In case of doped semiconduc- 
tors the example of Bi is of interest; interfaces in Bi- 
bicrystals of inclination type show superconductivit y up 
to 21 K, although Bi bulk is not a superconductor [20j . 
These two independently obtained indications, the pos- 
sible existence of high-temperature superconductivity in 
graphite / graphene and the special role of interfaces pit 
221 ] stimulated us to pursue the study of the trans- 
port properties of a bundle of internal interfaces in bulk 
HOPG samples. 

Transmission electron microscope (TEM) studies on 
HOPG samples revealed single crystalline regions of 



Bernal graphite of thickness in the c— axis direction be- 
tween 30 nm and ~ 150 nm[2l|. As example, we show in 
Fig. QJb) a typical TEM picture of the HOPG samples 
studied in this work. The different gray colors shown in 
Fig. QJb) indicate sightly different angle misalignments 
about the c-axis between each other and the existence 
of very well defined two-dimensional interfaces between 
them. We note that rotations up to 30° between the 
graphene layers from neighboring graphite regions has 
been seen by high resolution TEM in few layers graphene 
sheets [23j. The electrical response of those interfaces is 
the aim of our experimental study. 

Thin TEM lamellae have been prepared from a HOPG 
(ZYA grade, 0.3° rocking curve) bulk sample using a dual 
beam microscope (FEI Nanolab XT200). To avoid con- 
tamination and structural disorder the lamellae were cut 
after depositing a protective layer of 300 nm of tung- 
sten carbide using the electron gun (EBID) on top of 
the HOPG surface. The Ga + ion beam was used to cut 
the lamellae of thickness between ~ 300 and 800 nm 
in the a-b graphene plane direction and lengths up to 
~ 17 /im. After transferring the lamellae to an insulat- 
ing Si/SiN substrate, electron beam lithography followed 
by thermal evaporation of Pt/Au was used to make the 
four electrical contacts allowing us the measurement of 
the voltage drop of several interfaces in parallel to the 
graphene planes along a length ~ 2 . . . ~ 8 /im. More 
than ten lamellae were studied. Although we found qual- 
itatively similar behavior, the main difference in their 
characteristic transport properties depends on the used 
thickness, i.e. the observed Josephson-like I — V char- 
acteristics and the granular superconductivity behavior 
are observed at lower temperatures or even vanish the 
smaller the thickness (in the a, b plane) of the lamella. 
In this paper we present and discuss the results of four 
of them (L1-L4) with thickness (~ 500nm, ~ 800nm, 
~ 300nm, ~ 800nm), respectively. 

Four probe electrode configuration as shown in 
Fig.[lja), and also the van der Pauw configuration (con- 
tacts at lamella edges) were used to measure the tern- 
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FIG. 1. (a) Scanning Electron Microscopy (SEM) image of 
sample L3 on a Si/SiN substrate where the yellowish colored 
areas are the electrodes. A four points configuration has been 
prepared with the outer electrodes used to apply current and 
the inner ones to measure the voltage drop. As shown in the 
picture, the c-axis runs parallel to the substrate surface and 
normal to the current direction, (b) Transmission Electron 
Microscopy (TEM) image of one HOPG lamella. A different 
brightness corresponds to different orientation within the a — b 
plane of the crystalline regions of Bernal type with thickness 
between 30 and 200 nm. 
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FIG. 2. Temperature dependence of the voltage for four sam- 
ples measured with small input currents. A clear drop in the 
measured voltage is observed at 15 K< T < 150 K upon sam- 
ple. For the sample L4 the region near the onset of the voltage 
decrease is shown (second right y-axis). 



perature (T) dependence of the voltage (V) at constant 
input current (J) and the I — V characteristic curves. 
Figure [2] shows the T-dependence of the measured volt- 
age at I = 1 nA (L1,L3) and 100 nA (L2.L4). A clear 
drop in the measured voltage is observed, upon sample 
at different "critical" temperatures T c between ~ 15 K 
to ~ 150 K. This T c reflects the temperature above which 
the Josephson coupling between superconducting patches 
at some of the interfaces vanishes. At low enough temper- 
atures and currents zero resistance states were reached 
for LI, L3 and L4. Negative saturation voltages (L2) 
instead of zero are obtained for the van der Pauw con- 
figuration, which can be quantitatively explained using 
a simple Wheatstone bridge circuit and assuming two 
Josephson junctions with different T c 's as explained be- 
low. 

The observed T— dependence in Fig. [2] is determined 
mainly by the contribution of the internal interfaces, 
since the graphene layers within the crystalline regions 
show a semiconducting narrow-band- gap behavior [2J|. If 
the voltage drop is related to some kind of granular super- 
conductivity we expect that its temperature dependence 
is sensitively influenced by the applied current. Figures [3] 
(a) to (c) show clearly that the higher the input cur- 
rent the lower the transition temperature, revealing a 
semiconducting-like behavior at higher / and T . The ob- 
served behavior is compatible with the existence of gran- 
ular superconductivity (see, e.g., Ref. |25[ ) embedded in 
some of the internal interfaces. 

Further support to this claim is obtained from the 
I — V characteristics curves. Samples LI (Fig.^d)) and 
L3 (Fig. [3]Jf)) show typical Josephson behavior. In the 



case of L2 (Fig. E^e)) the curves were obtained with the 
van der Pauw configuration and a current and voltage 
paths that catch the answer of more than one Joseph- 
son junction in that sample. Taking into account the 
used current distribution in this case and assuming two 
Josephson junctions with different T c 's the rather exotic 
I — V curves can be quantitatively understood using the 
model of Ambegaokar and Halperin [26[ where the influ- 
ence of thermal fluctuations on the dc Josephson effect in 
a junction of small capacitance are taking into account, 
see Fig. SJa) . The same model can be successfully used 
to fit the I — V characteristics for sample L3 assuming 
only one Josephson junction, seeQJb). In the case of 
sample LI, a sharp jump in the current appears at the 
corresponding critical Josephson current I c where also a 
small hysteresis is observed, see Fig. [3](d). 

Considering the I — V results from all samples and 
configurations we obtain I C (T) shown in Fig. [5] in nor- 
malized units. The overall behavior is compatible with 
the temperature dependence expected for Josephson- 
j unctions where the normal barrier is given by ballistic 
graphene 27], the continuous line in Fig. O One may 
ask whether a Josephson coupling is possible through 
graphene layers and at large distances. Indeed, the work 
in Ref. 28 showed experimentally that the Josephson ef- 
fect is possible between superconducting electrodes sep- 
arated by hundreds of nanometers long graphene path. 
Therefore, the assumption of the existence of supercon- 
ducting regions Josephson coupled through graphene-like 
semiconducting paths at the observed interfaces appears 
reasonable. Because these superconducting regions at 
the interfaces are not homogeneously distributed in the 
used samples, upon the junctions' distribution, a noisy 
or "jumpy" behavior of the voltage is expected at low 
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FIG. 3. (a-c): Voltage vs. temperature at different input currents for samples L1-L3. The shadow region in (a) depicts the 300 
nV noise region of sample LI at low temperatures. This noise is intrinsic of the sample and vanishes after the application of a 
magnetic field, (d-f): I — V characteristics curves at different temperatures for the same samples. 
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FIG. 4. Current- Voltage characteristics at different tempera- 
tures for sample L2 (a) and L3 (b), this last in reduced coordi- 
nates where R is the normal state resistance. The continuos 
curves are fits to the model proposed in Ref. [2^ with the 
Josephson critical current I C (T) as the only free parameter. 



enough temperatures and currents due to phase and cur- 
rent path fluctuations. We have observed this behavior 
in samples LI, see Fig.[3£a), and L2, see inset in Fig.lBJa). 
This noisy behavior is intrinsic of the samples and can 
be suppressed by applying a magnetic field normal to the 
graphene planes, see inset in Fig. Eta), or increasing /. 
Future work should reveal the frequency spectrum of the 
noise and its relation to the Josephson junction arrange- 




Reduced Temperature (T/T c ) 



FIG. 5. Normalized Josephson critical current 7 C /7 C (0) vs. 
normalized T/T c obtained for different lamellae (each symbol 
corresponds to different a sample and/or configuration in the 
measurements). Critical currents between 55 nA and 5.5 fiA 
and critical temperatures from 15 K up to 175 K have been 
used. The continuous line is the theoretical curve taken from 
[2?| without any free parameter, assuming a short junction 
length. 



ment in the samples. 

An applied magnetic field is expected to be detrimental 
to the superconducting state. This effect can be due to 
an orbital depairing effect, i.e. a critical increase of the 
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FIG. 6. (a) Current- Voltage curves for sample L4 at 50K with 
and without magnetic field perpendicular to the graphene 
planes up to 8T. The inset shows the time dependence of 
the voltage in sample L2 for 250 nA input current at a con- 
stant temperature of 4 K, with no applied magnetic field 
(blue squares) and at 0.1 T field applied perpendicular to the 
graphene planes (red dots), (b) Current- Voltage characteris- 
tic curves for sample L3 at 2 K with magnetic field applied 
parallel to the graphene planes up to 8 T. The inset shows 
the temperature dependence of the voltage with 1 nA input 
current for sample L3 for magnetic fields from zero (blue cir- 
cles) up to 8 T (black squares) applied normal to the graphene 
planes. 



shielding currents, or at much higher fields due to the 
alignment of the electron spins, in case of singlet cou- 
pling. However, the possible effects of a magnetic field 
on the superconducting state of quasi two-dimensional 
superconductors or in case the coupling does not corre- 
spond to a singlet state are not that clear. For exam- 
ple, recently published experimental results [29j| in two 
different two-dimensional superconductors, including one 
produced at the interfaces between non superconduct- 
ing regions, show that superconductivity can be even en- 
hanced by a parallel magnetic field. In case the pair- 
ing is p— type[9( the influence of a magnetic field is ex- 
pected to be qualitatively different from the conventional 
behavior (30i l3ll] with even an enhancement of the super- 
conducting state at intermediate fields, in case the orbital 
diamagnetism can be neglected or for parallel field con- 



figuration. On the other hand, even for applied fields 
normal to the planes we expect much less influence of 
the orbital effect in case the London penetration depth is 
much larger than the size of the superconducting regions 
at the interfaces of our lamellae. If the superconduct- 
ing coherence length is of the order or larger than the 
thickness of the lamella then we expect a superconduct- 
ing (granular) behavior at lower T but a magnetic field 
may be less detrimental. We have studied the effects 
of magnetic fields applied parallel and normal to the in- 
terfaces on the transport characteristics in all measured 
lamellae. Upon sample the observed effects are from an 
usual detrimental, no effect at all or in some cases even 
a partial recovery of the superconducting state. 

Figure [5] shows the I — V characteristics of the sam- 
ples L4 (a) and L3 (b). For the relatively thick lamella 
L4 the magnetic field of 1 T applied normal to the in- 
terface planes vanishes the zero resistance state observed 
at zero field and at 50 K. At higher fields, however, the 
I — V curves show a recovery to the zero resistance state. 
A field applied parallel to the interface planes does not 
affect the curves. Recent studies on possible supercon- 
ductivity triggered by a large enough electric field ap- 
plied on multigraphene samples revealed also a kind of 
reentrance above a certain magnetic field applied normal 
to the graphene layers 32j. For the thinnest sample L3, 
however, a magnetic field applied in both directions has 
little or no effect on the I — V characteristics or V(T) up 
to 8 T, see Fig. EJb) and its inset. 

In conclusion, the transport characteristics of thin 
lamellae with tens of 2D interfaces between crystalline 
graphite regions reveal Josephson-like behavior with zero 
resistance states at low enough temperatures and input 
currents. Our results finally clarify that the origin for 
the metallic-like behavior as well as the giant magnetic 
field induced metal-insulator (MIT) transition measured 
in HOPG in the past 33|-35[ is related to the supercon- 
ducting properties of internal interfaces and it is not in- 
trinsic of the graphite structure. The existence of very 
high temperature superconductivity embedded at these 
interfaces is supported by recently done magnetization 
measurements [3fJ . 
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